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ABSTRACT 

Hinode discovered a beautiful giant jet with both cool and hot components 
at the solar limb on 2007 February 9. Simultaneous observations by the Hinode 
SOT, XRT, and TRACE 195 A satellites revealed that hot (~5x 10 6 K) and 
cool (~ 10 4 K) jets were located side by side and that the hot jet preceded the 
associated cool jet (~ 1 — 2 min.). A current-sheet-like structure was seen in 
optical (Ca u H), EUV (195 A), and soft X-ray emissions, suggesting that mag- 
netic reconnection is occurring in the transition region or upper chromosphere. 
Alfven waves were also observed with Hinode SOT. These propagated along 
the jet at velocities of ~ 200 km s~ x with amplitudes (transverse velocity) of 
~5-15 km s _1 and a period of ~ 200s. We performed two-dimensional MHD 
simulation of the jets on the basis of the emerging flux - reconnection model, 
by extending Yokoyama and Shibata's model. We extended the model with a 
more realistic initial condition (~ 10 6 K corona) and compared our model with 
multi-wavelength observations. The improvement of the coronal temperature 
and density in the simulation model allowed for the first time the reproduction 
of the structure and evolution of both the cool and hot jets quantitatively, 
supporting the magnetic reconnection model. The generation and the propaga- 
tion of Alfven waves are also reproduced self-consistently in the simulation model. 
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Introduction 



The solar ch romosphere has been known to be very dynamic fe.g.. iBrav fc Loughhead 



1974J ; IZirin Ill988l ). However, recent Hinode obse rvations (IKosugi et a 



is even more dynam ic than previously thought (jShibata et al.l 



2007 



2007) revealed that it 



Katsukawa et al. 2007 



Cirtain et al.l 120071 ) (see also special issue, Initial Results from Hinode, 2007, PASJ, Vol. 59, 
No. SP3). It h as long been observed that Ha jets c alled surges often occur in the chromo- 



sphere (see e.g. lRustlll967l; iKurokawa fc 



magn etic reconnection (IHeyvaerts et al. 



<Cawai 



19771 : 



1992) . Surges are believed t o be produced by 



Yokoyama fc Shibatalll995l . Il996l ; llsobe et al. 



20051 ) . which is an energy conversion mechanism from magnetic energy into thermal and ki- 
netic energies of plasma when two anti-parallel magnetic fields encounter and reconnect with 
each other. New chromospheric obser yations with the cal cium u H-line broad band filter 
of the Solar Optical Telescope (SOT; iTsuneta et al.l 120081 ) on board Hinode revealed that 
jets are ubiq uitous in the chromosphere and some of th e jets show evidence of magnetic 
reconnection (jShibata et al.l 120071 ; iKatsukawa et al.l 120071 ) . 



There have been s everal multi-wavelen gth observations of jets (e.g. X- r ay and Ha 



(ISchmieder et all Il995l : ICanfield et al.lll996h. Ha and E UV dChae et all Il999t iJiang et al. 
2007 ). and EUV and X-ray (j Alexander fc Fletcher 1^9991 )). These have shown that hot jets 



(X-ray or EUV jet) and cool jets (Ha surge or dark EUV jet) are almost co-spatial and co- 
temporal, and that they are dynamically connected to each other. Yokoyama and Shibata 
(1995, 1996) explained such a relation between hot and cool jets qualitatively by performing 
a resistive MHD simulation of an emerging flux model. However, their model was not fully 
realistic because of various computational difficulties; it was not possible to reproduce both 
emerging flux and jets self-consistently in one model, so they assumed an unrealistic initial 
temperature (~ 2.5 x 10 5 K) and density (10 12 cm -3 corona) for the corona. Because of this, 
the coronal Alfven speed in their model was lower than the actual value, so the jet velocity 
was also lower than the observed velocity of the Ha surge/jet. In this Letter, we succeed 
in modeling for the first time both the emerging flux and jet self-consistently with realistic 
coronal temperatures (~ 10 6 K) and densities and attempt to compare simulation results 
and observations quantitatively. 



2. The 2007 February 9 Giant Ca Jet Event 

The solar jet studied here occurred on the west limb of the Sun around NOAA active 
region 10940 on 2007 February 9 at 13:20 UT, which was associated with a weak soft X-ray 
brightening observed in GOES. The maximum height of the jet was ~ 14,000 km and its 
width ~ 6, 000 km. Figure la-lc shows Ca u H-line broad-band filter snapshot images taken 
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with SOT on board Hinode, which provide a high spatial resolution (0.2" or 150 km on the 
solar surface) and stable observation of the photosphere and chromosphere. The cold plasma 
jet, which we call the Ca jet, began to be ejected upward with inclination angle of 45 degree 
at 13:18 UT. The jet grew up to a cusp- or in verted- Y shaped stru cture, whose morphology 
was similar to that of coronal anemone jets (IShibata et al.lll994j ). The maximum upward 
velocity occurred 10 minutes later at 13:30 UT (also see Fig.2a). The jet was accelerated to 
~ 100 km s" 1 along the field line. The cool plasma subsequently fell down with a decelerated 
free-fall motion. 

Yohkoh observations have led to the interpretation of the anemone-shaped structure 
as a result of magnetic rec onnection between an emerging magnetic bipol e and a preexist- 
ing coronal uniform field (IShibata et al.l Il992l . 1 1994 : IShimojo et al.l Il996l ) . This interpre- 
tation has been suppor t ed by magnetohydrodynamic (MHD) simulations o f emerging flux 
(IHeyvaerts et al. 1 119771 ; lYokoyama fc Shibatal Il995l . 1 19961 : llsobe et al.l 120051 ) . By extending 
previous simulations the new simulations consider a more realistic initial condition, which 
allow a comparison with observations. Our initial condition of the simulations is ba- 
sically similar to that of Yokoyama & Shibata (1995), with the only difference being in 
coronal temperature and density; they were assumed to be 2.5 x 10 5 K and 10 12 cm -3 in the 



Yokoyama-Shibata model, but 10 6 K and 10 
parameters are much more realistic than theirs 



10 cm 3 



in our model. Note that our coronal 



As a result of our new two-dimensional MHD simulation, we found that a simulated 
jet is amazingly similar to the observed jet. The comparison among them is shown in 
Figure 1 (d-f are simulation results). These simulations are performed by solving the two- 



1 Initial coronal and chromospheric/photospheric temperatures are 100 and unity in unit of 10 K, respec- 
tively, and plasma (3 (= ratio of gas pressure to magnetic pressure) is 4 in the horizontal flux sheet just below 
the photosphere, 0.01-0.06 in the inner corona, and varies from 4 to 0.06 in the chromosphere/photosphere. 
Note that this simulation is based on an idealized model, in which there is no difference between chromo- 
sphere a nd photosphere in temperature. Numerical computation was carried out using the CIP-MOCCT 
scheme ( Kudoh et al. 19991 ) with total grid points of (2000 x 1000), whose grid size is .2 (<~ 40 km) and 
uniform in the whole computational box. We assumed anomalous resistivity ( Ugai 1985h . whose functional 
form is ij = for Vd < v c and rj oc (vd/v c — l) 2 for Vd > v c , where Vd = J/p, is the drift velocity, p is the 
mass density, J is the current density and v c is the threshold above which anomalous resistivity sets in. 
Typical magnetic Reynolds number VaL/t] at the location of reconnection is ~ 1000, where Va is the Alfven 
speed just outside the current sheet, L is the length of the current sheet, and rj is the magnetic diffusivity. 
The magnetic Reynolds number used in our simulation model is much lower than the actual number for 
solar corona, and the assumed functional form for the anomalous resistivity has not yet been physically 
established. However, it should be noted that this type of anomalous resistivity re produce well the fast 
recon nection model (e.g., Ugai 1985h and explosive energy release in solar flares (e.g.. lYokovama fc Shibata 



19941 ). 
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dimensional resistive MHD equations with uniform gravitational field and without thermal 
conduction and radiative cooling effects. We note that a jet is magnetically driven, not 
gas pressure driven, so the dynamics of the jet cannot be influenced by thermal conduction 
and radiative cooling effects even if the temperature is not fully realistic. We also note 
that evaporation effects due t o thermal conduction are not included in our model , so the 
evaporation-driven X-ray jets (IShimojo et al. Il200ll ; iMiyagoshi fc Yokoyama 1 120041 ) cannot 
be modeled in the Letter. As an initial state, we set hydrostatic plasma in the corona 
(f0 6 K), chromosphere/photosphere (f0 4 K), and convection zone with a horizontal flux sheet 
just below the photosphere and a uniform oblique magnetic field outside the flux sheet in 
the whole region. This flux sheet is unstable for the magne tic buoyancy instability (the 
Parker instability; IShibata et al.lll989l ; iMatsumoto et al.lll993l ). so the perturbed flux sheet 
is excited and emerges into the corona. While we performed a two-dimensional simulation, 
the emerging proce ss of magnetic flux would actually become mu ch more complex in three 
dimensions (see e.g. lFan. Abbett fc Fisher 1120031 ; llsobe et al.ll2005l ). In two dimensions, once 
reconnection occurs between the emerging flux and preexisting field, the field lines with a 
polarity opposite to that of the ambient field become connected to the ambient polarity 
regions, forming an inverted-Y shape or anemone jet. Because the time and spatial scales 
are arranged almost in the same way as in Figure la-c, the simulation model turns out 
to be able to explain observational facts very well. Reconnection creates multiple islands 
which confine cool, dense, chromospheric plasma in the current sheet, which are ejected at 
the Alfven speed va ~ 10 2 km s -1 (B/20G)(n/10 12 cm~ 3 ) _1 / 2 ~50-150 km s _1 if we assume 
coronal field B~20 G and density n~ f n ~ 12 cm -3 in the upper chromosphere, in agreement 
with the observations (shown in Fig.2a). These facts support a magnetic reconnection model. 
Note that, if reconnection occurs below the lower chrom osphere, shocks are f ormed in front 



of a jet due to rapid decrease of plasma density (e.g., IShibata et al.l Il982l ; iTarbell et al. 



1999|), which may eventually accelerate jets along magnetic field lines. However, in such 
a case, only cool jets are formed. In our case, it seems that reconnection occurred in the 
transition region or upper chromosphere, because not only cool jets but also hot jets are 
observed simultaneously (see section 3). 

Figure 2b shows the distance-time diagram of the position of the Ca jet (as represented 
by Ca intensity distribution), which revealed the oscillation of the jet with amplitude of 5-15 
km s _1 and period of 200 s. From the oscillation pattern at three different heights, we find 
that the oscillation propagates along the jet at 100-200 km s _1 . Since the jet is believed to be 
along the magnetic field, the propagation of the oscillation is the evidence of the propagating 
Alfven waves. The simulation model reproduced also the generation and propagation of 
Alfven waves (see plus symbols in Fig. 2b). Hinode observations revealed the existen ce of 
Alfven waves or Alfven-like oscillation (jOkamoto et al. 1120071 ; iDe Pontieu et al. 1120071 ). but 
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propagating Alfven waves have not been observed until now. Hence our observations are the 
first observational evidence suggesting the presence of propagating Alfven waves. 



3. Comparison between Hot/Cool components of the Jet and MHD simulation 



Our simulation shows that both hot and cool jets can be accelerated simultaneously 
by magnetic reconnection driven by emerging flux, as in Yokoyama & Shibata (1995). The 
plasma in the corona is heated to temperatures from a few million K to about 10 million 
K. This hot plasma can be observed as microflares and soft X-ray jets. At the same time, 
cool jets are accelerated if the reconnection occurs in the transition region or in the upper 
chromosphere, where cool plasma is situated near the reconnection point. According to the 
model, X-ray or EUV jets are seen as hot jets that are accelerated by the magnetic tension 
force of the reconnected field lines, and Hce surges or Ca jets are seen as cool jets that 
are accelerated by the slingshot effect due to the reconnection, which produces a whiplike 
motion. Furthermore, if a large amount of energy is injected into the upper chromosphere, 
the cold plasma above the energy injection point is heated up and evolves into evaporation 



( IShimojo et al. I l200ll ). although we cannot reproduce it because of the lack of thermal 



conduction in the energy equation. Such a coexistence of hot and cool jets is indeed confirmed 



by observations with t 



(Schmieder et al. 


1995; 


Canfield et al. 


1996 


Ko et al. 


20051 


Jiane et al. 


2007|). 



r e comparisons between associated X-ray or EUV jets and Ha surges 



1999 



Figure 3a-o shows the time evolution images of the jet with muft iwavelength observations 



Handy et al. 



ub et al. 2007) aboard Hinode, and the 



1999J) 195A filter. An X- 



taken with the SOT and X-Ray Telescope (XRT; iGo 
Transition Region and Coronal Explorer (TRACE; 
ray image before ejection shows a loop structure over the limb, at the footpoint of which a 
jetlike feature appeared at 13:09 UT, preceding other wavelengths. X-ray emission gradually 
increased its intensity and, just after 13:16 UT, an X-ray jet started to be ejected, and 
cusp- or inverted-Y shaped structure was formed. An EUV jet was ejected from a looplike 
bright emission patch at the moment when the patch reached its maximum intensity. The 
cold plasma ejection, such as the Ca jet and a dark EUV jet in absorption (~ 10 4 K), was 
delayed by 1-2 minutes relative to the X-ray and EUV brightening. The maximum upward 
velocity of the cool jets occurred about 10 minutes later than the X-ray spike (also see 
Fig. 2a). As a result of spatial co-alignment, the EUV jet was identified with the X-ray jet, 
whereas the dark EUV jet appeared to be a counterpart to the Ca jet. The X-ray jet and 
the Ca jet were ejected side by side with each other, which is the same feature as shown 
in the simulation results of Figure li. These observational facts indicate the following three 
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features: (1) The X-ray jet (5 x 10 6 K) and the EUV jet (~ 10 6 K) are likely to have the same 
physical origin, and the Ca jet (10 4 K) and the dark EUV jets (unknown temperature, ~ 10 4 
K) also have the same origin. (2) The X-ray and the Ca jets are different kinds of plasma 
ejection, implying dynamically connected hot and cool plasma ejections along different field 
lines. (3) The X-ray jet precedes the Ca jet by 1-2 minutes compared with their maximum 
intensity time, and some hot structure was also observed with X-ray emission 10 minutes 
earlier than jet ejection. 

We note that the delay of the cool jet may not be due to a cooling effect, because the 
delay time is much shorter than the cooling time. Furthermore our simulation without a 
cooling term reproduces both hot and cool jet very well. According to the simulation, the 
current sheet is formed between the corona (hot and thin plasma) and the emerging flux with 
chromospheric de nsity and temperature (coo l and dense plasma). Hence the reconnection is 



very asymmetric (IPetschek fc Thorne 1 119671 ). In such a case, the low-density part becomes 
a hot and fast jet, while the high-density part forms a cool and slow jet, because the local 
Alfven speed is high (low) in the low (high) density part. This is why a hot jet preceded a 
cool jet; i.e., the hot jet reaches higher altitude earlier than the cool jet. 



4. Summary and Discussion 



MHD simulation results (Fig.ld-i) reproduce remarkably well the dynamics and struc- 
ture of cool and hot jets and their relative timings. The only structure that the simulations 
cannot explain is the existence of an X-ray bright point (probably an unresolved loop) in 
Fig.3m-o. However, MHD simulations in Figure 1 are two-dimensional, and hence discrep- 
ancy between observations and simulations may be explained by three-dimensional effects . 



In fact, the reconnec tion model (jYokoyama fc Shibatalll995l . Il996l ; IShibata et al.lll992l . 11994 



Shimojo et al.lll996l ) predicts the formation of not only jets but also loops in a separated 
place, which can explain man y X-ray observations showing that bright poin ts (loops) are 
situated separately from jets (IShibata et al.lll992l . Il994 ; IShimojo et al.lll996l ). In our case, 
such bright loops may be situated just in front of jets because of three-dimensional projection 
effects. 

On the other hand, detailed comparison between simulations and observations suggests 
that the current sheet structure in Figure le may be visible as one of legs in the inverted-Y 
shaped structure. It is interesting to see that EUV and X-ray loops seem to be situated 
along the same leg, i.e., possibly corresponding to the current sheet. 



We note here that an X-ray brightening at the footpoint cannot be seen in EUV emission 



- 7- 



but in absorption. This may be interpreted as either the three-dimensional effect that the 
EUV emission from the X-ray source is covered by cool plasma ejection or the temperature 
effect such that the temperature in the X-ray source is too high to emit enough EUV emission. 
(Note that usually X-ray loops cannot be seen in EUV images.) 

It is also noted that both observations and simulations show the generation of Alfven 
waves when magnetic reconnection occurs: the jet undergoes apparent motion perpendicular 
to the jet direction at amplitudes of 5-15 km s" 1 (see Fig. 2b). This velocity is also com- 



parable to that observed f or polar X-ray jets (ICirtain et al.l 12007). Th ese Alfven waves are 



generated by reconnection (jYokoyamalll998l ; iTakeuchi &: Shibata 



20011 ) , and may contribute 



to the heating and accele ration of solar wind when the magnetic field is open (IParker 111988 
Suzuki fc Inutsuka 1 120051 ) . 
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Fig. 1. — Comparison between the Ca jet and simulated jets based on the reconnection 
model, (a)-(c) Ca II H broadband filter images of the Ca jet on 9 February 2007, taken with 
Hinode/SOT. (d)-(f) Two-dimensional distributions of density (logio P\ color map) in units 
of 10~ 7 g cm -3 in the simulated jets, whose times correspond to those of the observed jets 
in the top panels. The length is in units of 200 km and the times are in units of 20 s. (g-i) 
Two- dimensional distribution in density (logio P\ color maps) in units of 10 -7 g cm -3 and 
temperature (logio T; color map) in units of 10 4 K of magnetic fields (B; lines) and velocity 
vectors (v; arrows) in the simulated jets. 
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Fig. 2. — (a) Height-time diagram of intensity (running differences) along the Ca jet taken 
with Ca n H broadband filter of Hinode SOT. The numbers above the curves denote velocity 
(in km s _1 ). The plus symbols show the height of the simulated jet shown in Fig. 1. The 
simulation reproduced well the observed motion of the jet. (b) Distance-time diagram of Ca 
intensity across the jet, revealing the oscillation of the jet at amplitude of 5-15 km s -1 and 
period of 200 s, at three locations (heights). It is found that the oscillation propagates from 
lower altitude at slit 1 position to higher altitude at slit 2 and 3 positions. Since the jet 
is considered to be parallel to magnetic field lines, this is the evidence of the Alfven wave 
propagating at the velocity 100-200 km s -1 . The simulated jet also shows the generation 
and propagation of Alfven waves along the cool jet, which is indicated by the plus symbols 
(the position of the jet). 




Fig. 3. — Comparison of multi- wavelength observations, (a-e, f-j, k-o) Time evolution images 
of the Ca jet taken with Hinode/SOT Ca II H-line broad band filter, the EUV jet with 
TRACE/1951 filter and the X-ray jet with Hinode/XRT Al poly filter. 



